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A  Global  Thermospheric  Density  Model 
Based  on  Satellite  Accelerometer  Data 


1.  INTRODUCTION 

Extensive  measurements  of  the  lower  thermosphere  neutral  density  have  been 
derived  from  the  accelerometers  on  the  Atmosphere  Explorer  -C,  -D,  and  -E, 
and  S3-1  satellites.  An  empirical  Accelerometer  Density  Model  (ADM)  has  been 
developed  from  these  data.  This  model  permits  organization  of  the  large  den¬ 
sity  data  base  in  a  form  suitable  for  geophysical  interpretation  and  for  comparison 

12  3 

with  other  data.  As  in  previous  studies,  ’  '  the  least-squares  method  of  multi¬ 
ple  linear  regression  is  used  to  extract  and  analyze  the  density  data.  Density  vari¬ 
ations  in  the  altitude  region  140  to  240  km  are  represented  as  a  function  of  geomag¬ 
netic  activity,  latitude,  local  time,  day  of  year  and  solar  flux  (60  £  F  £  130  units). 
The  major  advantage  of  the  present  model  is  improved  formulations  for  the  latitude 
variations  in  geomagnetic  activity,  local  time  and  semiannual  effects. 


(Received  for  publication  13  January  1982) 
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Knowledge  of  atmospheric  density  variability  is  required  for  satisfactory  de¬ 
sign,  operation,  and  orbit  prediction  of  low -altitude  AF  mission  satellites.  Exist¬ 
ing  models  of  the  lower  thermosphere,  particularly  below  250  km,  are  deficient 
since  they  are  primarily  based  on  data  at  higher  altitudes.  Physical  and  mathe¬ 
matical  difficulties  limit  extrapolation  of  these  data  to  lower  altitudes;  namely, 
the  assumption  of  diffusive  equilibrium  is  not  always  valid,  and  the  lower  thermo¬ 
sphere  boundary  conditions  are  poorly  known  and  variable.  This  report  compares 

"  ‘  4  5 

the  ADM  results  with  the  three  most  commonly  used  models:  J71,  J77,  and 

£ 

MSIS.  Results  of  this  study  provide  a  quantitative  estimate  of  improvements  re¬ 
quired  in  properties  of  current  thermospheric  models  in  order  to  account  for  the 

observed  variations.  Additional  low -satellite  altitude  accelerometer  data  obtained 
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during  moderate  and  high  solar  flux  conditions  are  becoming  available.  ’  These 
results  will  be  utilized  in  order  to  extend  the  present  formulations  to  all  solar  flux 
conditions. 


2.  DATA  SELECTION  AND  COVERAGE 


Table  1  summarizes  orbital  characteristics  and  dates  of  acquisition  for  each 
q 

satellite.  Marcos  et  al  have  provided  a  complete  description  of  the  experiments 
and  the  data  base.  Over  three  million  density  measurements  were  obtained.  This 
amount  is  too  large  to  handle  practically  and  much  larger  th  n  required  to  define 
the  parameters  to  be  fitted.  Therefore,  a  data  base  was  constructed  by  selecting 
values  at  5-km  intervals  from  140  to  240  km.  At  each  5-km  interval,  the  density 
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Table  1.  Satellite  Orbital  Characteristics  and  Data  Acquisition  Periods 


Satellite 

Launch 

_  - 

End  Elliptical 
Data  Acquisition 

Inclination 

(degrees) 

Initial 

Perigee 

(km) 

Initial 

Apogee 

(km) 

AE-C 

Dec  73 

Nov  74 

68° 

156 

4000 

S3  - 1 

Oct  74 

May  7  5 

97° 

160 

4000 

AE-D 

Oct  75 

Jan  7  6 

90° 

156 

3800 

AE-E 

Nov  75 

Nov  76 

2  0° 

157 

3000 

value  was  selected  using  a  second-order  fit  to  three  points:  the  one  closest  to  the 
5-km  altitude  step  and  the  two  adjacent  points.  As  a  result,  data  from  the  four 
satellites  are  represented  by  148,  000  density  values. 

Typical  results  from  the  data  base,  namely,  individual  density  values  at 
180  km,  are  shown  versus  day  of  year  in  the  second  from  the  top  frame  of  Fig¬ 
ure  la,  b,  and  c.  The  10.7-cm  solar  flux,  ratio  to  the  J7i  model  and  the  3-lir 
planetary  index  are  also  given.  The  1974  data  (Figure  la)  are  from  AE-C 
(1  January  to  30  October),  and  from  AE-C  and  S3-1  (3l  October  to  3i  December). 
The  1975  data  (Figure  lb)  are  from  S3-1  to  15  May;  from  AE-D  (13  October  to 
6  December),  and  from  AE-D  and  AE-E  (7  December  to  31  December).  The  1970 
data  (Figure  lc)  are  from  AE-D  and  AE-E  to  30  January,  and  from  AE-E  alone 
to  15  November.  Distribution  of  the  complete  data  base  as  a  function  of  day  of 
year,  altitude,  latitude,  local  time,  solar  flux  and  geomagnetic  activity  has  been 

Q 

given  by  Marcos  et  al  for  each  individual  satellite  data  set.  They  are  summa¬ 
rized  and  discussed  briefly  here,  with  ail  four  sets  of  data  combined. 

The  distribution  of  data  points  with  altitude  is  shown  in  Figure  2a.  Data  below 
150  km  were  obtained  only  during  the  AE-C  and  AE-E  perigee -lowering  maneuvers. 
The  latitude  distribution  (Figure  2b)  peaks  between  ±  20°  (41  percent  of  data  base) 
because  of  the  contribution  of  the  AE-E  data.  However,  there  are  56,  000  points 
from  20°  to  60°  I  and  31,  000  points  from  1 60°  to  90°  1 .  These  are  adequate  to  per¬ 
mit  analysis  of  latitudinal  variations  in  density.  Excellent  local  time  coverage 
was  provided  by  AE-C  and  -E,  whereas  AE-D  and  S3-1  data  were  mainly  confined 
to  local  times  near  2200  hours  and  0800  hours  (Figure  2c).  Data  were  obtained 
during  the  declining  portion  of  solar  cycle  20  and  the  beginning  of  solar  cycle  2  1. 
Figure  2d  shows  that  flux  values  greater  than  100  were  exceeded  only  4  percent  of 
the  time.  A  wide  range  of  geomagnetic  conditions  was  encountered  (see  Fig¬ 
ure  2e).  Geomagnetically  active  conditions  (assuming  Kp  >  3)  occur  for  40  per¬ 
cent  of  the  data.  Latitudinal  coverage  as  a  function  of  day  of  year  is  shown  in 
Figure  3.  Position  of  the  upleg  and  downleg  latitude  for  data  at  the  typical  altitude 
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lb.  Accelerometer  Data  Base  at  180  km  for  1975.  Latitude,  density 
),  ratio  to  J7 1  model  and  K_  as  a  function  of  time 
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Figure  lc.  Accelerometer  Data  Base  at  180  km  for  1976.  Latitude,  density 
(g/cm^),  ratio  to  J7  1  model  and  K  as  a  function  of  time 
&  P 
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Figure  2a.  Histogram  of  Data  Distribution  as 
a  Function  of  Altitude 
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Figure  2b.  Histogram  of  Data  Distribution  as 
a  Function  of  Geographic  Latitude 
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Figure  2c.  Histogram  of  Data  Distribution  as 
a  Function  of  Local  Time 
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SOLAR  FLUX 


Figure  2d.  Histogram  of  Data  Distribution  as 
a  Function  of  Solar  Flux 
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Figure  2e.  Histogram  of  Data  Distribution  as 
a  Function  of  Geomagnetic  Activity 
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of  200  km  1b  shown  for  each  satellite.  Coverage  from  AE-E  Is  evident  as  the  data 
points  cycling  between  ±  20°.  Atmosphere  Explorer -D  provides  measurements 
extending  to  high  latitudes  for  northern  hemisphere  winter  and  southern  hemisphere 
summer  conditions;  S3-1  follows  a  similar  track  but  continues  to  northern  hemi¬ 
sphere  spring  conditions.  The  AE-C  data  include  the  high  latitude  northern  hemi¬ 
sphere  spring  and  southern  hemisphere  fall  conditions. 

One  advantage  of  this  data  set  is  that  all  measurements  were  obtained  by  the 
same  technique  and  have  essentially  the  same  accuracy  (see  Section  4.  5).  Hence 
no  one  satellite  data  set  will  have  undue  weighting  influence  on  the  fit  of  specific 
parameters. 

3.  FORMULATION  OF  THE  ACCELEROMETER  DENSITY 

MODEL (ADM) 

A  multiple  linear  regression  model  is  used  to  approximate  the  dependent  vari¬ 
able,  atmospheric  density  by  a  linear  combination  of  independent  variables.  The 
independent  variables  considered  in  this  study  (at  each  5 -km  interval)  are  the  vari¬ 
ations  associated  with  solar  flux,  geomagnetic  activity,  day  of  year,  local  time, 
and  latitude.  Terms  were  chosen  to  represent  known  atmospheric  variations  in 
a  relatively  simple  manner.  All  latitude  dependencies  are  described  in  geographic 
coordinates.  Test  runs  indicated  that  the  fit  to  the  data  was  not  sufficiently  im¬ 
proved  to  warrant  the  complexity  of  a  combination  of  geographic  and  geomagnetic 
coordinates. 

Table  2  contains  a  list  of  the  independent  variables.  Variations  related  to 
solar  activity  are  approximated  by  the  81 -day  centered  average  of  the  observed  ' 
10.  7-cm  solar  radio  flux  (Fjq  7)  and  by  the  difference  Fj0  7  -  Fj0  7  using  a  one  day 
lag  for  the  daily  F^Q  7  value.  Geomagnetic  activity  variations  are  represented, 
respectively,  by  the  average  of  the  four  3-hr  K  values  (K^j.,1  and  the  eight  K 
values  (Kp24^  that  occurred  before  the  time  of  the  density  measurements,  by  a 
latitude  dependence  for  K  hv  K  with  a  6-tir  lag  (K  and  hv  a  latitude- 

pl2  p  p6  _0 

dependent  time  lag.  This  lag,  derived  hv  correlation  analvsis  of  data  in  la 

latitude  intervals  is  of  the  form: 

K  At  =  3  +  A-  |  |  0]  -  60 1  hours  (1) 

p  12 

where  9  =  geographic  latitude.  This  is  an  average  result.  The  response  time  is 
believed  overestimated,  since  data  are  not  available  on  a  continuous  basis,  or  even 
on  an  every-orbit  basis,  and  the  measurements  are  not  always  in  the  localized 
region  of  initial  heating.  The  latter  situation  will  particularly  overestimate  the 
high-latitude  time  delay.  However,  Eq.  (1)  is  consistent  with  the  measured  data. 


Table  2.  Terms  Used  in  Regression  Analysis 


1. 

F 

10.7 

2. 

F  -  TT* 

10.7  10.7 

3. 

RP12 

4. 

KP24 

5. 

KP6 

6. 

Cos  8  Kp 

7. 

Sin  0  *S>12 

8. 

K  At 

P 

9. 

COS  ut  +  dj 

10. 

cos  2ut  + 

11. 

cos  3ut  +  <^3 

12. 

cos  8  (cos  ut  +  d4) 

13. 

sin  8  (cos  ut  + 

14. 

cos  8  (cos  2wt  +  ^g) 

15. 

sin  8  (cos  2u  t  +  <t> 7> 

16. 

cos  8  (cos  3ut  +  <ba) 

O 

17. 

sin  8  (cos  3ut  +  dg) 

18. 

cos  *ss 

19. 

cos  at  + 

2  0. 

cos  2at  +  $  j  j 

21. 

sin  8  (cos  at  + 

22. 

sin  8  (cos  2at  +  ^ j3) 

23. 

te| 

Solar  Flux 


Geomagnetic 

Activity 


Local 

Time 


Day  of 
Year 


J  Latitude 


Variations  associated  with  day  of  year  are  represented  by  sine  and  cosine  functions 
having  periods  of  1  and  l/2  year  and  a  nonsymmetric  latituue  dependence.  Local 
time  effects  are  represented  by  terms  of  the  form  cos  (net)  +  d  where  <t>  -  phase  and 
n  1,2, 3  and  by  latitude-dependent  components  of  the  form  (cos  8)  cos  nut  and 
(sin  8)  cos  nut.  Also  a  cos  xgg  term  was  used  to  represent  a  subsolar  bulge  with: 


18 


cos  x  =  cos  6  cos  t  cos  6  +  sin  6  sin  6  (2) 

ss 

where 

X  -  solar  zenith  angle 
6  =  latitude 
t  =  local  time 
6  =  solar  declination  angle. 

A  pure  latitude  variation  was  included  by  the  absolute  value  of  the  satellite  geograph¬ 
ic  latitude. 


4.  RESULTS 

With  use  of  the  terms  in  Table  2,  the  stepwise  multiple  linear  regression 
analysis  was  carried  out  for  each  5-km  interval  data  set  between  140  and  240  km. 

As  in  previous  analyses,  for  example,  Forbes  and  Garrett,  the  process  of  adding 
terms  in  each  case  was  terminated  when  the  Student's  t-value  fell  below  the  95  per¬ 
cent  confidence  level.  Results  for  successive  applications  of  the  regression  analy¬ 
sis  to  data  at  200  km  are  given  in  Table  3.  This  table  lists  the  functions  (F.)  in 
order  of  entry  into  the  regression  analysis,  regression  coefficient  (a.),  their 
standard  deviations  (cr.l,  and  the  correlation  coefficients  (R)  for  the  final  regression 
equation.  Although  there  are  altitudinal  differences  in  the  relative  importance  and 
the  entry  of  some  of  the  coefficients,  the  200-km  data  were  selected  as  being  typi¬ 
cal.  Density  values  are  calculated  by: 

p  =  a  +  E.  a.  F.  (3) 

o  111 

4.1  Solar  Activity  Variation* 

The  major  energy  source  affecting  the  thermosphere  is  the  solar  ultraviolet 
flux  between  30  and  2400  A.  The  global  mean  heating  due  to  the  solar  EUV  heat 
source  is  about  4  X  10**  erg -sec  1  above  12  0  km  and  10 erg -sec  *  above  100  km. 
A  correlation  of  density  with  the  10.  7-cm  radio  flux  monitored  at  NRC  Ottawa. 
Canada  is  a  feature  of  all  current  models.  This  correlation  does  not  necessarily 
represent  the  complex  physical  mechanisms  of  the  interactions  between  the  solar 
EUV  flux  and  the  thermosphere.  However,  since  F  }  7  is  found  to  reflect  vari¬ 
ations  in  the  thermospheric  energy  input,  it  provides  a  readily  available  but  imper¬ 
fect  indicator  of  solar  activity. 

The  density  response  to  the  average  solar  flux  F  in  the  empirical  model  (ADM) 
is  examined  by  letting  F  =  F  (so  that  F-F  =  0).  For  a  unit  increase  in  F  the  data 
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Table  3.  Regression  Analysis  Results  at  200  km 


i 

Fi 

— 

a. 

i 

ai 

R 

1. 

V2 

1.  544  E-14 

9.  962  E-16 

0.  5449 

2. 

F10.  7 

2. 063  E - 15 

5.251  E- 17 

0.  647  3 

3. 

cos  (2a  +  109.  30) 

2.  045  E-14 

4. 569  E-16 

0.  7091 

4. 

sin  8  cos  (a  +  159.  34) 

2. 519  E-14 

7.210  E-16 

0.7376 

5. 

F10.  7  '  F10.  7 

9.  033  E- 16 

3.441  E-  17 

0.7548 

6. 

cos  8  cos  (ut  +  14.  94) 

1.  526  E-14 

6. 409  E-16 

0.  7720 

7. 

cos  (a  +  28.  2  5) 

1.  144  E-14 

5.  234  E-16 

0.  7876 

8. 

%24 

6.298  E-  15 

6.  407  E-16 

0.  7901 

9. 

sin  8  cos  (2ut  +  0.70) 

1.  958  E- 15 

8. 190  E-16 

0.  7  922 

10. 

cos  8  cos  (3ut  +  4.  17) 

8.  710  E-15 

1.  468  E-15 

0.  7930 

11. 

sin  8  cos  (2a  +  48.  10) 

1.  911  E-15 

8.217  E-16 

0.  7937 

12. 

cos  8  Kpi2 

-8.  095  E-15 

7. 970  E-16 

0.7  942 

13. 

1  8  | 

-2. 001  E - 16 

2.758  E - 17 

0.7968 

14. 

cos  Out  +  0.  38) 

4.706  E-15 

1.  173  E-15 

0.7972 

15. 

K  At 

P 

1.  187  E-15 

4.  654  E-16 

0.7977 

16. 

sin  8  K 

-9.  652  E-16 

1.  958  E-16 

0.7978 

p  12 

17. 

cos  (2ut  +  1.28) 

1. 859  E-14 

1.  504  E-15 

0.  7980 

18. 

cos  8  cos  Out  +  7.  19) 

2. 090  E-14 

1.761  E-15 

0.  8012 

19. 

sin  8  cos  (ut  +  0.  8 1) 

3. 422  E-15 

8.  185  E-16 

0.  8014 

20. 

cos  (ut  +  10.  54) 

1. 588  E-15 

4.636  E-16 

0.  8017 

0  =  geographic  latitude 
u  =  2  ir/24 
a  =  2  ir/  3  6  5 .  2 

local  time  phases  are  in  hours 

day  of  year  phases  are  in  days 

a  =  1.722  1  -14 
o 


analysis  Indicates  a  density  increase  varying  from  0.25  percent  at  140  km  to 
1.  50  percent  at  240  km.  These  computations  are  made  for  the  March  equinox 
(day  81)  at  the  equator  with  local  time  =  1500  hours  and  »  0.  These  default  val¬ 
ues  are  used  in  the  remainder  of  this  section.  The  ADM  results  are  represented  by 
the  solid  straight  line  in  Figure  4;  also  shown  are  the  J71,  J77,  and  MS1S  increases 
per  unit  F  increase  computed  at  20-km  intervals.  At  140  km  the  MS1S  value,  not 
shown  on  the  plot,  is  -0.  04  percent  (a  density  decrease  is  predicted).  At  all  alti¬ 
tudes  the  ADM  shows  the  largest  response  to  F.  The  J77  model  response  is  about 
0.  1  percent  lower,  showing  the  closest  agreement.  The  actual  mass  density  values 
as  a  function  of  F  are  shown  for  the  altitude  of  200  km  in  Figure  5.  Note  that  the 
ADM  has  the  lowest  density  value  at  F  =  60.  For  the  range  F  =  60  to  130  units, 
the  average  density  increases  per  unit  F  are  1.  08  percent  (ADM),  0.  98  percent 
(J77),  0.33  percent  (J7l)  and  0.72  percent  (MSIS). 


140  160  180  200  220  240 


ALTITUDE  (Km) 

Figure  4.  Percentage  Density  Increase  per  Unit  Increase  in  Average  Solar  Flux 
(F)  Plotted  vs  Altitude  (Day  *  81,  I  T  =  1500  hr,  Latitude  -  0,  K  =0,  F  =  751 

Density  variability  at  200  km  with  respect  to  the  daily  solar  flux,  F,  is  ex¬ 
amined  in  Figure  6.  For  this  case,  F  varies  from  60  to  130  units  whereas  F  is 
held  constant  at  75  units.  (Again  the  ADM  has  the  lowest  density  values  for  low 
values  of  F.  )  The  ADM  increases  by  0.45  percent  per  F  unit.  The  other  three 
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AVERAGE  SOLAR  FLUX  (F) 


Figure  5.  Density  Variations  at  200  km  as  a  Function  of  the  Mean  Solar 
Decimetric  Flux  F.  for  Day  *  81,  LT  =  1R00  hr,  Latitude  =  0°  and  K  -  o 

P 


Figure  6.  Density  Variations  at  200  km  as  a  Function  of  F-F;  for  F  =  7i 
Day  *  81,  I  T  -  lnOO  hr.  Latitude  =  0°  and  K  =0 


models  give  the  values  0.  57  percent  (J77),  0.  32  percent  (J71),  and  0.  18  percent 
(MSIS).  For  this  case,  the  ADM  data  predict  a  variation  that  is  about  halfway 
between  the  J77  and  J71  model  values. 

It  is  pointed  out  that  these  results  are  obtained  during  the  declining  portion  of 
solar  cycle  20  and  the  beginning  of  solar  cycle  21.  Linear  extrapolation  of  the 

present  model  to  higher  solar  flux  conditions  may  not  be  valid.  This  can  be  ex- 

7  8  10 

amined  with  ROCA  and  SETA  data.  Also,  Hinteregger  has  shown  that  solar 

EUV  radiation  correlates  more  closely  with  7  during  solar  cycle  21  than  it 

did  for  solar  cycle  20.  Therefore,  in  any  model,  caution  is  required  in  basing 

density  predictions  on  results  derived  from  a  different  solar  cycle.  This  emphasizes 

the  need  for  a  solar  activity  index  based  on  direct  determination  of  the  solar  EUV 

heating. 

4.2  Geomagnetic  Variation* 

Rapid  global  fluctuations  in  density  with  a  typical  duration  of  one  to  several 
days  are  connected  with  geomagnetic  variations.  The  energy  input  results  from 
solar  wind/earth's  magnetic  field  interactions  and  involves  Joule  heating  and  parti¬ 
cle  precipitation  along  magnetic  field  lines.  During  geomagnetically  quiet  conditions, 

17  - 1 

the  energy  input  due  to  these  high  latitude  processes  is  of  the  order  5  X  10  erg  -  sec 
This  energy  input  can  increase  by  a  factor  of  10  or  more  during  large  geomagnetic 
storms.  The  planetary  3 -hr  K  index  used  as  an  indicator  of  geomagnetic  activity 
does  not  necessarily  represent  physical  mechanisms  responsible  for  density  vari¬ 
ations.  It  is  derived  from  a  network  of  poorly  distributed  magnetometer  stations 
(ten  between  43  and  60°N  and  two  between  37  and  43°S).  Hence  it  is  particularly 
limited  when  describing  data  at  high  or  low  latitudes.  As  with  Fjq  however, 

Kp  is  presently  the  only  index  available  for  systematic  geomagnetic  activity 
correlations. 

All  models  show  similar  trends  that  are  associated  with  geomagnetic  storms. 

The  amplitude  of  the  density  increase  generally  depends  on  the  strength  of  the  geo¬ 
magnetic  storm  as  indicated  by  the  index.  The  atmospheric  response  occurs 
first  at  auroral  zone  latitudes  and  then  at  low  latitudes,  with  a  delay  of  typically 
6  to  8  hours.  This  implies  that  the  energy  transferred  to  the  thermosphere  through 
auroral  processes  at  high  latitudes  is  distributed  later  through  dynamic  processes. 
Recent  composition  measurements  and  theoretical  studies  have  shown  that  atmos¬ 
pheric  winds  driven  by  the  high-latitude  heat  source  are  involved  in  global  redis¬ 
tribution  of  energy  and  mass.  This  wind-induced  dij  usion  at  high  latitudes  results 
in  a  decrease  in  the  light  atmospheric  constituents  (O  and  He),  whereas  the  heavy 

10.  Hinteregger,  H.  E.  (1979)  Development  of  solar  cycle  21  observed  in  EUV 
spec*— um  and  atmospheric  absorption,  J.  Geophys.  Res.  84:1933. 


constituents  (N2  and  A)  increase  in  response  to  the  temperature  increase.11  Con¬ 
sequently,  the  response  of  neutral  density  is  expected  to  be  latitude  and  composition 
dependent. 

The  density  variations  derived  for  the  geomagnetic  effect  exhibit  an  altitude 

dependence.  Figure  7  compares  the  increase  in  density  per  unit  rise  in  K  as  a 

P 

function  of  altitude  predicted  by  the  ADM  and  J7l,  J77,  and  MSIS  models.  Per¬ 
centages  are  computed  as  1/ 6  of  the  increase  in  density  going  from  =  0  to 
Kp  *  6.  The  data  are  for  F  =  F  =  75,  day  =  81,  latitude  =  0  and  LT  =  1500  hours. 
Regression  analysis  revealed  increases  from  4.4  percent  at  140  km  to  P.  7  percent 
at  240  km  for  the  ADM.  On  the  whole,  J7 1  gives  closest  agreement  to  this  result. 
Both  J7  1  and  MSIS  show  less  variability  than  ADM.  However,  J77  predicts  a  con¬ 
stant  value  of  11.7  percent  at  all  altitudes.  This  is  a  consequence  of  partitioning 
the  geomagnetic  effect  into  three  parts,  with  the  "density  wave"  affecting  the 
equatorial  density.  At  higher  latitudes,  the  J77  model  predicts  an  altitude -dependent 


ALTITUDE  (km) 


Figure  7.  Percentage  Density  Increase  per  Unit  K 
Increase  as  a  Function  of  Altitude 


11. 


Mayr,  H.G.,  Harris.  I. 
atmosphere  dynamics. 


,  and  Spencer,  N.W.  (1P78)  Some  properties  of  upper 
Rev.  Geophys  Space  Phvs.  l(i:53P. 
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Figure  8  shows  the  increase  in  density  at  200-km  altitude  as  a  function  of  K^. 
At  values  of  >  6,  the  ADM  values  are  shown  as  a  dashed  line.  Extrapolation  of 
a  linear  model  is  not  valid  at  the  high  values  where  the  geomagnetic  effect  be¬ 
comes  highly  nonlinear.  (The  low  density  value  at  *  0  in  the  MSIS  data  is  due  to 
a  continuity  problem  in  the  model  rather  than  an  assumed  nonlinear  effect. )  Limi¬ 
tation  of  the  ADM  results  to  K  =6  only  excludes  about  1.  5  percent  of  the  K  values 

P  P  12 

experienced,  based  on  the  present  data  set  and  long-term  occurrence  probabilities. 


Figure  8.  Density  Variation  at  200  km  as  a  Function 
of  Kp  (Day  =  81,  I.T  =  1500  hr.  Latitude  =  0,  F  -  !•'  -  75) 


A  significant  latitude  dependence  was  derived  in  the  regression  analysis.  Fig¬ 
ure  9  compares  the  percentage  increase  in  density  per  unit  versus  geographic 
latitude  predicted  by  the  various  models.  Data  are  presented  for  the  altitude  of 
2  00  km.  At  the  equator,  the  ADM  enhancement  is  7. 3  percent;  it  increases  to 
13.  1  percent  at  -90°  and  to  11.4  percent  at  +90°.  The  small  pole-pole  difference 
is  within  the  standard  deviation  of  each  enhancement  and  cannot  be  unambiguously 
related  to  a  hemispheric  asymmetry  in  the  geomagnetic  effect.  The  shape  of  the 
ADM  curve  and  its  range  of  values  agrees  best  with  MSIS.  Since  higher  order  terms 

12.  Cage,  A.  L.  and  Zawalick,  E.J.  (1972)  A  Discussion  of  Geomagnetic  Indices 
Kp  and  Ap,  1932  to  1971,  AF<  RL-72-0693.  AD  756828. 
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J77 
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MSIS 


Figure  9.  Percentage  Density  Increase  per  Unit  Kp 
as  a  Function  of  Geographic  Latitude  1 


were  not  analyzed  in  the  ADM  formulation,  it  is  not  possible  to  evaluate  the  compli¬ 
cated  J77  latitudinal  response.  Averaged  over  all  latitudes,  however,  J77  agrees 
well  with  ADM. 

At  all  altitudes,  the  term  with  the  highest  correlation  coefficient  was 
Table  3  shows  that  at  200  km  this  term  accounted  for  29.  7  percent  of  the  variation 
in  the  data.  Regression  analysis  tends  to  smooth  the  short-term  or  transient  vari¬ 
ations  in  the  data;  examination  of  Figure  1  shows  individual  magnetic  storm  effects 
with  response  related  to  a  specific  jump  in  Kp.  Statistical  analyses  (unpublished) 
of  the  present  accelerometer  data  also  indicate  geomagnetic  variations  that  are  a 

function  of  local  time,  season,  and  hemisphere.  These  will  be  examined  in  a  fur- 

7  8 

ther  study  incorporating  the  more  recent  ROCA  and  SETA  accelerometer  results. 
A  major  limitation  of  all  present  models  is  their  inability  to  predict  accurately  the 
phase  and  amplitude  of  geomagnetic  variations.  However,  this  difficulty  is  being 
alleviated  by  new  analyses  in  both  the  MSIS^’  and  J77^  models. 


References  13  to  15  are  listed  on  page  42. 


2fi 


4.3  Local  Time  Variation! 


Density  structure  as  a  function  of  local  time  has  been  extracted  in  previous 
studies  of  accelerometer  density  data.  Marcos  et  al9  statistically  studied  AE-E 
data  at  160,  200,  and  240  km,  and  found  characteristics  of  atmospheric  tides.  At 
160  km,  two  nearly  equal  maxima  were  observed  near  0730  and  2300  hours  re¬ 
spectively;  at  240  km,  a  broad  maximum  was  centered  near  1400  hours.  Forbes 
16 

and  Marcos  extended  this  AE-E  analysis.  They  determined  that  the  transition 

from  predominantly  semidiurnal  to  predominantly  diurnal  tidal  effects  occurs  at 

180  ±5  km  and  that  a  substantial  summer/winter  difference,  consistent  with  the 
17 

Forbes  model,  occurs  in  the  phase  structure  of  the  semidiurnal  component  of 

1 8 

density.  Using  data  from  all  four  satellites,  Forbes  and  Marcos  showed  that 
the  phase  of  the  diurnal  variation  varies  with  latitude  shifting  to  later  local  times 
from  the  equator  to  midlatitudes. 

The  ADM  density  structure  as  a  function  of  local  time  is  shown  in  Figures  10a, 
b,  and  c  for  the  altitudes  160,  200,  and  240  km,  respectively.  For  each  case,  den¬ 
sity  values  are  calculated  assuming  day  =  81,  latitude  =  0,  F  =  F  =  75  and  Kp  -  0. 

In  Figure  10a  the  ADM  results  show  the  dominance  of  the  semidiurnal  tide.  This  re¬ 
sult  is  in  remarkable  agreement  with  MSIS.  Table  4  shows  that  the  phases  of  the 
maxima  and  minima  for  these  models  are  within  l/2  hour  for  Uiree  cases  and  with¬ 
in  one  hour  for  the  fourth  case.  Ratios  of  the  morning  and  evening  maxima  to  the 
early  morning  and  afternoon  minima  are  compared  in  Table  5.  The  ADM  ampli¬ 
tudes  are  in  the  range  14.  5  ±  1.  5  percent,  whereas  those  of  MSIS  are  within  15.  0 
±1.5  percent.  Both  J7l  and  J77  show  only  a  diurnal  variation  with  their  maximum 
(amplitude  of  about  29  percent)  near  1500  hours.  This  time  corresponds  to  near- 
minimum  densities  in  ADM  and  MSIS.  At  200  km  (Figure  10b),  ADM  predicts  a 
broad  maximum  between  1100  and  2000  hours  as  the  diurnal  variation  becomes  more 
prominent.  This  result  agrees  well  with  the  MSIS  local  time  behavior,  but  not  with 
the  well-defined  diurnal  bulges  of  the  J7  1  and  J 77  models.  For  each  model  the  ratio 
of  maximum  to  minimum  density  is  as  follows:  1.20  (ADM),  1.  15  (MSIS),  1.62  (J71), 
and  1.  56  (J77).  All  models  show  the  dominance  of  the  diurnal  tide  at  240  km 
(Figure  10c).  A  phase  of  1330  hours  for  the  maximum  is  obtained  for  ADM.  The 
J7  l  and  J77  phases  are  at  1400  and  1530  hours,  respectively,  whereas  the  MSIS 
data  maximize  at  1600  hours.  The  amplitude  of  the  ADM  variation  (1.  54)  is  between 
thai  of  the  MSIS  (1.26)  and  the  J7  1  ( 1.  97)  and  J77  (1.  95)  models.  Based  on  the 

16.  Forbes,  J.  M.  ,  and  Marcos,  F.  A.  (1979)  Tidal  variations  in  total  mass  density 

as  derived  from  the  AE-E  Mesa  experiment,  J.  Geophys.  Res.  84:31. 

17.  Forbes,  ■!.  M.  (1978)  Tidal  variations  in  O,  09,  N?,  Ar,  He  and  H,  J.  Geophys. 

Res.  83:3691. 
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Figure  10a.  Density  Variation  at  160  km  as  a  Func¬ 
tion  of  Local  Time  (Day  *  81,  Latitude  =  0,  F  =  P  =  75, 
Kp  =  0) 
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Figure  10b.  Density  Variation  at  200  km  as  a  Func¬ 
tion  of  Local  Time 
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Figure  10c.  Density  Variation  at  240  km  as  a  Func¬ 
tion  of  Local  Time 


Table  4.  Phase  of  Semidiurnal  Variation  at  100  km 


Maxima  Minima 

(hr)  (hr) 

ADM  8  1/2  21  l/2  2  l/2  15  i/2 

MSIS  9  21  3  1/2  lo 


w 
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Table  5.  Amplitude  of  Semidiurnal  Variation  at  100  km 


Morning  Max 

Evening  Max 

r  .... 

Morning  Max 

Evening  Max 

Afternoon  Min 

Afternoon  Min 

Morning  Min 

Morning  Min 

ADM 

1.  140 

1.  130 

1.  lr.o 

1.  150 

MSIS 

1.  135 

1.  157 

1.  143 

_ 

1.  104 

_ 

present  analysis,  the  J 71  and  J77  models  significantly  overestimate  the  amplitude 
of  the  diurnal  bulge  at  the  altitudes  studied.  The  cos  \  term  (Eq.  2)  representing 
the  subsolar  bulge  did  not  enter  the  regression  analysis  at  any  altitude  between  140 
to  240  km.  Examination  of  the  ADM  density  distribution  as  a  function  of  latitude 
(for  LT  =  1500  hr)  during  the  solstices  shows  maximum  density  values  at  auroral, 
rather  than  subsolar,  latitudes  in  the  summer  he  mi:  phere.  This  is  assumed  to 
be  a  consequence  of  high  latitude  heat  sources. 

4.4  Annua)  Variations 

Density  variations  with  annual  and  subannual  periodicities  as  distinguished 
from  periodic  changes  connected  with  the  rotation  of  the  earth  and  sun  are  a  regular 
feature  of  the  upper  atmosphere.  They  are  referred  to  as  the  "semiannual  vari¬ 
ation"  and  are  generally  characterized  in  current  models  by  equinoctial  maxima 
and  solstitial  minima.  At  the  present  time  the  cause  of  the  semiannual  variation 

i  q 

is  not  completely  understood.  Ching  and  Chiu  1  have  shown  that  EUV  heating  has 

proper  characteristics  to  be  important  at  altitudes  above  150  km.  Mayr  and 
20 

Volland  have  shown  that  sources  in  addition  io  EUV  heating  are  required  to  ac¬ 
count  fully  for  this  variation.  High  latitude  heat  sources,  including  auroral  heating 

’>  1 

related  to  the  semiannual  component  of  the  occurrence  of  magnetic  storms,”  have 
been  invoked.  More  data  and  theoretical  analysis  are  required  to  develop  a  com¬ 
prehensive  model  of  the  semiannual  variation. 

The  semiannual  density  variation  at  200  km  extracted  from  the  accelerometer 
data  at  the  equator  and  north  and  south  poles  is  shown  in  Figure  11a.  Kesults  are 
for  F  *F  =  75,  LT  =  1500  and  Kp  =  0.  Data  at  corresponding  latitudes  are  given  for 
MSIS  and  J77  models  in  Figures  lib  and  11c,  respectively.  Data  for  .'7]  are  not 

19.  Ching,  B.K.,  and  Chiu,  Y.T.  (1972)  Annual  and  sub-annual  effects  of  EUV 

heating-I.  Harmonic  analysis,  Planet.  Space  Sci.  20:1745. 

20.  Mayr,  H.G.,  rnd  Volland,  H.  (1972)  Theoretical  model  for  the  latitude  depend¬ 

ence  of  the  annual  and  semiannual  variations,  J.  C.eophys.  Kes.  77  :0774. 

21.  Russell,  C.  T.  ,  and  McPherron,  R.  L.  (1973)  Semiannual  variation  of  geomag¬ 

netic  activity,  J.  Geophys.  Res.  78:92. 
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DAY  COUNT 

Figure  lie.  Density  at  200  km  vs  Day  of  Year  for  MSIS 

shown,  but  are  essentially  the  same  as  for  J77.  Phases  at  the  three  latitudes  are 
summarized  in  Table  6.  Although  there  are  variations  in  the  month  of  occur¬ 
rence  of  some  maxima  and  minima,  the  columns  are  labeled  "Jan  min.  ,  "  "Apr 
max.  ,  "  "  Julv  min. "  and  "Oct  max.  "  for  convenience.  In  each  case,  the  phases 
and  their  variations  with  latitude  are  in  very  good  agreement.  Although  the  J77 
January  minimum  occurs  20  days  before  that  of  ADM  and  MSIS  at  -90°,  the  curves 
are  quite  flat  for  this  case.  Therefore,  the  discrepancy  is  not  significant.  A  most 
striking  feature  of  the  amplitude  of  the  ADM  semiannual  variation  is  the  latitudinal 
behavior  of  the  October  maximum.  In  the  southern  hemisphere,  this  maximum  is 
larger  than  the  April  maximum.  At  the  equator,  the  two  maxima  are  nearly  equal. 
At  the  north  pole,  the  October  maximum  amplitude  is  very  small.  The  details  of 
this  behavior  are  qualitatively  similar  for  the  other  two  models  depicted  in  Fig¬ 
ures  lib  and  11c.  The  difference  between  the  hemispheres  occurs  because  the 
seasonal  effect  is  in  phase  with  the  annual  effect  in  the  southern  hemisphere  and  in 
opposite  phase  in  the  northern  hemisphere.  To  analyze  further  the  derived  ampli¬ 
tude  structure,  ratios  have  been  calculated  at  -90°,  0°,  and  +90°  for  each  model 
for  Tour  cases:  ( 1  >  Apr  max/July  min,  (2)  Oct  max/July  min,  (3)  Apr  max/Jan 
min  and  (4)  Oct  max;Jan  min.  These  results  are  summarized  in  Table  7.  At  the 
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Table  6.  Phase  of  Semiannual  Variation  at  200  km 


Latitude 

Model 

Jan  Min 
(day) 

Apr  Max 
(day) 

July  Min 
(day) 

Oct  Max 
(day) 

ADM 

40 

90 

195 

310 

-90° 

J77 

20 

90 

2  00 

310 

MSIS 

40 

80 

190 

ADM 

20 

100 

190 

mm 

0° 

J77 

15 

100 

2  05 

300 

_ _ 

MSIS 

20 

100 

190 

290 

ADM 

5 

110 

2  05 

270 

+90° 

J77 

10 

100 

210 

290 

M^IS 

10 

11 0 

190 

280 

- 

Table  7.  Amplitude  of  Semiannual  Variation  at  2  00  km 


Latitude 

Model 

Oct  Max 
July  Min 

_ 

Apr  Max 
Jan  Min 

Oct  Max 
Jan  Min 

ADM 

_ 

1.  03 

1.  20 

-911° 

J77 

1. 09 

1.  19 

MSIS 

1.  01 

1.  19 

ADM 

1.  36 

1.  33 

1.  18 

1.  16 

0° 

'77 

1.29 

1.  34 

1.  17 

1.22 

MSIS 

1.24 

1.31 

1.20 

1.27 

ADM 

1.  33 

1.  10 

1.  50 

1.23 

+90° 

J77 

1.25 

1.23 

1.28 

1.25 

MSIS 


1.29 


I.  35 


1.  58 


1.  36 


equator,  all  three  model  amplitudes  are  in  excellent  agreement  with  values 
within  ±1.  5  percent  for  cases  2  and  3,  and  within  ±6  percent  for  cases  1  and  4. 

At  -90°,  the  ADM  amplitude  structure  agrees  more  closely  with  MSIS.  For  the 
north  pole  data,  the  comparison  results  are  less  consistent.  All  models  show  an 
amplitude  of  29  ±4  percent  for  case  1.  The  ADM  amplitude  agrees  better  with 
MSIS  for  case  3  and  with  J77  for  case  4.  The  ADM  Oct  max/July  min  (case  2) 
amplitude  (10  percent)  is  lower  than  both  J77  (23  percent)  and  MSIS  (35  percent). 
Therefore,  the  ADM  October  maximum  exhibits  a  strong  hemispheric  dependence 
with  a  significant  decrease  in  amplitude  from  the  south  to  north  pole. 

Global  density  distributions  at  160,  200,  and  240  km  are  summarized  as  a 
function  of  latitude  and  day  of  year  in  Figures  12a,  b,  c,  respectively.  These 
contour  maps  are  for  the  conditions  F  =  F  =  75,  LT  =  1500  and  Kp  =  0.  Numbers 
on  each  plot  indicate  a  density  value  relative  to  the  average  density  (for  all  condi¬ 
tions)  at  that  altitude.  Examination  of  these  figures  shows  that  the  characteristics 
of  the  semiannual  variation,  particularly  the  low  amplitude  of  the  northern  hemi¬ 
sphere  October  maximum,  are  similar  at  all  three  altitudes. 


Figure  12a.  Contour  Map  of  Density  at  160  km  vs  Latitude  and  Day  of  Year,  for 

F  *  F  =  75,  I  T  -  1500  hr  and  K  -  0 

P 
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Composition  data  at  200  km  from  MSIS  at  -90°,  0  and  +90°  are  shown  in  Fig¬ 
ure  13.  Nitrogen  CNj)  and  atomic  oxygen  (O)  display  different  behavior  as  a  func¬ 
tion  of  latitude  and  day  of  year.  While  atomic  oxygen  has  a  pronounced  semiannual 
variation  at  +90°,  Nj  has  one  maximum  in  May.  However,  composition  effects  are 
not  distinguishable  in  accelerometer  data.  Composition  measurements  from  mass 
spectrometers  on  AE-D  should  aid  in  explaining  the  relative  concentrations  of  O 
and  N2  leading  to  the  low  total  density  values  for  the  northern  hemisphere  fall  maxi¬ 
mum.  However,  the  cause  of  this  phenomenon  and  whether  it  is  peculiar  to  low 
solar  flux  conditions  require  additional  data  and  theoretical  analysis. 
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Figure  13.  Composition  Data  at  200  km  vs  Day  o  1  Year  for  MSIS 


4.5  Statistical  Evaluation  of  Models 

To  investigate  the  over-all  accuracy  of  models,  statistical  comparisons  to  the 
data  base  were  made.  The  first  and  second  moments  of  R  (ratio  of  measured  data 
to  model)  are  determined  respectively  by: 

*  •  I  -if  «> 

i=l 
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and 


a 


(Bi  -  R)2 

N-l 


1 


(5) 


where  R  is  the  mean  value,  R^  is  the  ratio  of  the  i^  density  measurement  to  the 
model,  N  the  number  of  points  and  a  the  standard  deviation. 

The  measured  density  values  from  each  individual  set  of  satellite  data  (at  all 
altitudes)  were  first  statistically  analyzed  in  relation  to  the  J71,  J77,  and  MSIS 
models.  Table  8  summarizes  values  of  N,  R  and  a  for  each  satellite.  The  mean 
values  are  within  ±2  percent  for  each  set  with  respect  to  the  J71  and  J77  models 
and  within  ±3  percent  for  the  MSIS  model.  This  demonstrates  the  consistency  of 
the  data  sets  for  modeling  purposes.  Further,  it  implies  that  the  average  drag 
coefficient  for  S3-1  did  not  differ  significantly  from  that  of  the  AE  satellites;  S3-1 
was  trapezoidal -shaped  while  the  AE's  were  squat  cylinders.  For  all  four  satel¬ 
lites  a  constant  value  of  Cj-j  =  2.2  was  used.  Standard  deviations  in  Table  8  fall  in 
the  range  0.  14  to  0.  18  for  all  cases.  The  MSIS  provides,  for  the  average  conditions 
of  this  data  base,  a  slightly  better  depiction  of  atmospheric  variability.  The  per¬ 
centage  standard  deviation  for  each  model,  averaged  over  the  four  sets  of  satellite 
data  are;  MSIS  =  14.7,  J71  =  15.  1  and  J77  =  15.8. 


Table  8.  Ratio  of  Total  Mass  Density  to  Model  Density  for  Satellites  S3-1,  AE-C, 
-D  and  -E  for  Altitudes  140-240  km 


N 

J71 

J77 

MSIS 

R 

0 

R 

0 

R 

0 

S3-1 

27,315 

1.  06 

0.  15 

1. 04 

0.  15 

1.  00 

0.  14 

AE-C 

57,716 

1.  09 

0.  16 

1.05 

0.  17 

1.  05 

0.  15 

AE-D 

28,668 

1.  05 

0.  16 

1.00 

0.  15 

0.  98 

0.  16 

AE  -E 

34,364 

1.  08 

0.  18 

1.  03 

0.  18 

_ 

1.  01 

0.  15 

Evaluation  of  R  and  a  for  ADM  was  made  at  each  5 -km  interval,  using  the  com¬ 
bined  four  data  sets.  Table  9  shows  a  comparison  of  ADM  and  the  other  three 
models  at  the  altitude  of  20 0  km.  The  ADM  percentage  standard  deviation  is 
1.9  percent  less  than  that  of  MSIS,  2.  7  percent  less  than  J71.  and  3.  4  percent  less 
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Table  9.  Comparison  of  Model  Density 
Ratios  for  All  Satellite  Data  at  200  km 


R 

a 

ADM 

1.  00 

0.  12 

J71 

1.  09 

0.  16 

J77 

1.  04 

0.  16 

MS  IS 

1.01 

0.  14 

Total  number  of  points  =  8747 

than  J77.  The  calculated  standard  deviation  for  ADM  must  presently  be  viewed  as 
an  internal  test  of  the  consistency  between  the  mathematical  representation  and  the 
data  set  used  for  model  development.  The  test  of  the  model  is  its  ability  to  repre¬ 
sent  data  not  used  in  its  construction.  An  evaluation  of  a  preliminary  version  of 
this  model  with  available  data  obtained  over  a  20-day  period  from  a  high  inclination 
satellite  did  show  results  consistent  with  those  of  Table  9.  More  extensive  tests 
are  planned  using  a  larger  data  set. 

5.  CONCLUSIONS 

Density  variations  in  the  lower  thermosphere  related  to  solar  flux,  geomagnet¬ 
ic  activity,  local  time,  latitude  and  the  semiannual  effect  have  been  determined  from 
the  S3-1  and  AE-C,  -D,  and  -E  accelerometer  data  base,  using  multiple  linear 
regression  analysis.  The  resulting  Accelerometer  Density  Model  is  believed  to 
provide  a  better  representation  of  the  low -altitude  satellite  aerodynamic  environ¬ 
ment  than  other  global  density  models.  The  major  results  of  this  study  are: 

1.  The  density  is  more  variable  with  the  solar  cycle  dependent  EUV  flux  as 
represented  by  7  than  predicted  by  current  models. 

2.  Geomagnetic  activity  variations  exhibit  a  latitude  dependence  with  a  maxi¬ 
mum  at  high  latitudes. 

3.  The  subsolar  (diurnal  tide)  bulge  is  dominated  at  altitudes  below  ~2  00  km 
by  the  semidiurnal  tide.  Also,  density  values  maximize  in  the  summer  hemisphere 
at  auroral  latitudes  rather  than  at  low  latitudes  in  the  vicinity  of  the  subsolar  point. 
This  is  consistent  with  a  high-latitude  heat  source  that  is  effective  even  during  geo- 
magnetically  quiet  times. 
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4.  A  strong  hemispherical  asymmetry  for  the  October  maximum  character¬ 
izes  the  semiannual  variation.  The  October  maximum  has  its  greatest  amplitude 
at  the  south  pole  and  its  amplitude  decreases  from  the  south  pole  to  the  north  pole. 

5.  Statistically,  the  ADM  standard  deviation  is  about  2  to  3  percent  lower  than 
that  of  the  other  models  tested. 

The  accelerometer  data  base  has  been  used  to  develop  an  improved  model  of 
lower  thermospheric  density  variations.  This  model  is  applicable  to  low  solar  flux 
conditions  and  is  based  largely  on  data  at  low  and  middle  latitudes.  Data  obtained 
from  the  ROCA  and  SETA  accelerometer  experiments  on  high  inclination  satellites 
have  been  obtained  during  periods  of  moderate  to  high  solar  flux.  Incorporation  of 
these  results  into  the  present  data  base  will  improve  the  accuracy  and  extend  the 
applicability  of  the  present  model. 
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